Penetration of the magnetic field into the twinning plane in the type I and II 

superconductors 
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It is demonstrated that in the type I and II superconductors with low-transparent twinning 
planes (TP) the penetration of external parallel magnetic field into the region of the twinning plane 
can be energetically favorable. In the type I superconductors the twinning planes become similar 
to Josephson junctions and the magnetic field penetrates into the center of the TP in the form 
of soft Josephson-like vortices. This leads to increase in the critical magnetic field values. The 
corresponding phase diagram in the parameter plane "temperature - magnetic field" essentially 
differs from the one obtained without taking the finite value of the magnetic field near the TP into 
account. Comparison between obtained phase diagrams and experimental data for different type I 
superconductors can allow to estimate the value of the TP transparency, which is the only fitting 
parameter in our theory. 

PACS numbers: 74.25.Dw, 61.72.Mm, 74.50.+r, 74.78.Na 



The phenomenon of the twinning plane superconduc- 
tivity (TPS) has been the subject of intensive investiga- 
tions during last three decades (see Ref . [l| for review) . A 
twinning plane (TP) may produce more favorable condi- 
tions for the superconducting nucleation compared with a 
bulk crystal and a superconducting layer localized on the 
TP can appear even above the bulk critical temperature 
T c . Recently the interest to the physics of twins in super- 
conductors was renewed since it was shown experimen- 
tally that TP affect the properties of many relatively new 
superconductors which belong to the pnictide family 2 - - — . 
In these superconductors twinning planes can enhance 
locally the superfiuid density^ or influence the vortex 
pinning 7 -^. In particular it was demonstrated^ that in 
Ba(Fei_ x Co x )2As2 twinning planes repulse vortices and 
act as strong barriers for vortex motion. Thereupon the 
theoretical investigations of the magnetic properties of 
the TP in superconductors are of current importance. 

Twinning planes can effectively screen parallel mag- 
netic field which leads to the increase in the value of the 
critical field as a function of temperature. The corre- 
sponding phase diagrams H C (T) for the absolutely trans- 
parent TP were studied theoretically within the phe- 
nomcnological Ginzburg-Landau formalism both for the 
type and the type II superconductors^. In particu- 
lar for the ultra type I superconductors it was shownii - — 
that for small but finite values of the Ginzburg-Landau 
parameter k the magnetic field penetration into the su- 
perconducting area leads to the corrections to the TPS 
free energy proportional to k(" +1 )/ 2 (n = 0,1,...) and 
in practice only the term oc k 1 / 2 plays an important 
role while the terms of the order k 3 / 2 and higher can 
be neglected since they weakly contribute to the result- 
ing H C (T) diagram. The resulting dependencies H C (T) 
do not contain any fitting parameters, which allowed 
their quantitative experimental verification for concrete 
superconductors^^. In some case s 16 ' 17 the local en- 



hancement of superconductivity may occur near the sam- 
ple surface. The u ppe r critical field for this situation was 
considered in Ref. Il8l . However for the type I-Sn the su- 
perconductivity in magnetic field should appear as I or- 
der transition (except very narrow region near the critical 
temperature) and this field corresponds to the overcool- 
ing of the normal phase. 

Practically for TP with finite electron transparency the 
standard Ginzburg-Landau free energy functional should 
be generalized by inserting additional term which breaks 
the requirement of the order parameter continuity at 
TPi 19 i 20 The influence of the finite TP transparency on 
the upper critical magnetic field in the type II super- 
conductors was analyzed in Ref. [U At the same time 
for superconductors of the I type it was found that in 
case of low transparent TP the essentially asymmetric 
distributions of the order parameter relative to the TP 
can become energetically more favorable than symmet- 
ric onesi^ 2 - The striking prediction of this paper is that 
under certain conditions the order parameter is nonzero 
only at one side of the TP while at another side it should 
be zero. 

We would like to point out that all theoretical results 
which have been obtained up to now arc based on the 
assumption that superconducting nucleus localized near 
TP is screening the magnetic field effectively: it was be- 
lieved that in the type I superconductors the magnetic 
field can penetrate only into the region which is far from 
the TP (the magnetic field value at the TP is exponen- 
tially small) while in the type II superconductors the 
magnetic field value has its minimum at the TP. 

In the present paper we show that for both type I and II 
superconductors with low-transparent TP the magnetic 
field can fully penetrate into the center of the twinning 
plane and the corresponding state is energetically favor- 
able. For the type II superconductors this fact results in 
small corrections to the magnetic susceptibility only. At 
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the same time for the type I superconductors the mag- 
netic field penetration into the center of the TP leads 
also to essential changes in the dependence of the crit- 
ical magnetic field on temperature due to the negative 
contribution to the free energy, which has the order of k. 
Note that the obtained solutions have lower energy than 
the ones found in Ref. [22|. At the same time the corre- 
sponding profiles of the order parameter are symmetric 
relative to the TP. 



Let us consider a bulk superconducting sample with a 
single twinning plane at z = 0. The external magnetic 
field H is assumed to have only the y component. We 
choose the corresponding vector potential in the form 
A x (z) = Hz so that the order parameter ip depend only 
on x and z. We will use the standard Ginzburg-Landau 
free energy functional to describe the local enhancement 
of the superconductivity on the T P 11 ' 21 
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where a = a(T — T c ), B = rot A, p is a phenomeno- 
logical constant describing the finite transparency of the 
TP, V± = ip(x,y,±0) and the value £ s will be defined 
below. Let us also introduce the temperature dependent 
coherence length as £(T) = h / \/ '4ma(T - T c ) . 

For the further analysis it is convenient to rewrite the 
functional ([T]) in dimcnsionlcss variables 



1 = J - 1' Hs = Hc{t = 6 = e(< = 

~ ib x z 

ljj s = ip (t = -1), ijj = — , x = —, z = — , 

H s kH s ^ s 4 s 8tt 
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Here T s is the critical temperature of the superconductor 
with the twinning plane (T s > T c ), ipo(T) = y/\a\/b and 
H C (T) = Y / 47ra 2 /& are the values of the order parameter 
wave function and the critical filed for the bulk super- 
conductor without twinning plane; K = mcyb/y^Treh is 
the Ginzburg-Landau parameter. In what follows we will 
omit the tildes since we will consider only dimensionlcss 
expressions. Then the functional for the free-energy per 
unit length along the TP reads as 



G = G S dxdz < 2 
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Let us start with the case of the absolutely opaque TP 
(r = oo) in the ultra type I superconductor (k -C 1) and 
consider the order parameter wave function in the form 
ip(x,z) = ip(z)exp{iO{x)} . Varying the functional (O by 
ip* and A for z^O we obtain the equations 



d%tp + tip + tp 3 + ip(A/V2-d x 6) =0, 



d\A = (V2^ 2 /k 2 ) (a/V2 - d x e 



(4) 



with the boundary conditions d z tp+ = —tp+ and d z ip- 
<P-. Note that the system (j4]) has the first integral 



(d z ip) 2 -tip 2 
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The modulus of the order parameter on the TP is given 
byli <fl = (1 - t) + ^J{l-t) 2 -h 2 . In what follows 
we will assume that the magnetic field penetration into 
the center of the TP weakly affects the order parameter 
modulus near the TP. The validity of this assumption 
will be discussed below. Also we would like to mention 
that the order parameter decay length I is proportional 
to i -1 / 2 while the magnetic field penetration length is 
proportional to k/\/1 — t and for n <C 1 and (1— t) ^ k 2 
it is much less than I. This allows one to consider the 
constant order parameter value in the region near the 
center of the TP where the local magnetic field b = nd z A 
is nonzero. We will assume that b(z = ±oo) = h and 
b(z = 0) = ho, where ho < h. Then from the solution of 
the equation (jj) on the vector potential we obtain 



b(z) = /i exp — ! o x 6 = —j= 
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The resulting correction to the free energy value per 
unit length along the x axis has the form G = Gq + 
G s (2k/<po) (ho — 2hh ), where G is the free energy ob- 
tained without taking the penetration of the magnetic 
field to the TP into accounti±*H 
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From the obtained dependence G(h ) one can see that 
its minimum corresponds to the case ho = h and the 
resulting free-energy value is 



G = G -G s (2 K / lf0 )h 2 



(7) 
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FIG. 1: (Color online) Phase diagram of the type I supercon- 
ductors with absolutely opaque TP in the parallel magnetic 
field. Penetration of the magnetic field into the TP leads to 
increase in the critical magnetic field h c (t) (red solid curve) 
compared with the one calculated with the assumption of ex- 
ponentially small magnetic field on the TP^i (blue dashed 
curve). In the inset we show the fragment of these depen- 
dencies in more detail. In our calculations we took k — 0.13 
corresponded to tin. Also critical magnetic fields of the su- 
perconductor with k = (black solid curve) and for a bulk 
superconductor (black dashdot curve) are shown. 




FIG. 2: (Color online) Phase diagram of the type I super- 
conductors with low-transparent TP (r = 25) in the parallel 
magnetic field. Penetration of the magnetic field into the TP 
is energetically favorable for temperatures below the point 
where the dependence of the critical Josephson magnetic field 
h c j(t) (green dotted curve) crosses the critical field depen- 
dence calculated with the assumption of the exponentially 
small magnetic field on the TP^i (blue dashed curve) . In this 
temperature region the resulting dependence of the critical 
magnetic field is shown with red solid curve. In our calcula- 
tions we took k = 0.13 (tin). 



Thus it is energetically favorable for magnetic field to 
penetrate fully into the TP center. Then solving numer- 
ically the inequality G < we obtain the temperature 
dependence of the critical magnetic field h c (t), which is 
shown in Fig. 1 (red solid curve). 

Note that the magnetic field (|6|) does not affect the 
value of the order parameter at the TP. Indeed substi- 
tuting the corresponding vector potential at z = to the 
first integral and one can obtain that at the TP the 
order parameter ip 2 (z = 0) = 2(1 —t) — ft. 2 /Vo = fo- 

The expression for the local magnetic field ([BJ allows 
one to calculate the correction to the magnetic moment 
of the TP, which has the form (here we consider only the 
correction due to the magnetic field ([6]) and restore the 
dimension of the expression) 

oo 

AM, = l^ dz = pL. ( 8 ) 

(J) J 4-7T 2TTtp Q ^ ' 

— oo 

Comparing this value in the limit h — > with the 
diamagnctic moment Mmi^ due to the expulsion of 
the magnetic field from the superconducting region 
we obtain that the ratio ri(j\ = AMm/M^m = 

k (In k -1 ) -1 \A/ 2 (! _ *)• Practically r)^ < 1 for all 
temperatures where the TPS has the I type. Indeed, for 
example, for tin with n = 0.13 the condition rjtj-\ ~ 1 
gives (1 — t) < 2 • 10~ 3 . For such small values of 
(1 — t) the obtained results are not applicable since for 



(1 — t) < 2.7 • 10~ 2 the superconducting phase transition 
is of the II typei. 

Now let us generalize the obtained results for the case 
of TP with finite but small transparency (we will assume 
that r _1 <C 1). For finite r the TP is similar to the 
Josephson junction with corresponding Gibbs free en- 
ergy G r /G s oc r^ 1 + 0(r~ 2 ). Note that the boundary 
conditions for the order parameter tp at z = should 
be modified in order to describe the transparent TP— so 
that 

d z ip+ = -i/>+ + r~ 1 (ip+ - ip-), , s 

<9 Z V- =V- + r _1 (V+-^-)- [> 

Note that the corrections to the order parameter mod- 
ulus |^| due to the changes in the boundary conditions 
have the order of r _1 . The corresponding corrections to 
the free energy have at least the order of r~ 2 and can 
be neglected. Thus in what follows we will assume that 
for low-transparent TP the spatial distribution of the or- 
der parameter module (p(z) is the same as in case of the 
absolutely opaque TP. Then it is easy to obtain the ex- 
pression for the Josephson free energy of the TP, which 
has the form G r /G s = (4<^g/r) / (1 — cosA6>) dx, where 
A6 = 6+- 6_ and 9± = 9(x, ±0). 

It is naturally to expect that the magnetic field 
can penetrate into the TP center in the form of soft 
Josephson-like vortices. Indeed the magnetic field at the 
TP is defined by the phase difference on the two sides of 
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the TP and has the form ho(x) = (ipo/\/2)d x A8 while 
dimensionless Josephson current through the TP can be 
written as 



Jz{x) = ^M sin A9 



(10) 



(here we use the value cHg/A-K^s as the unit of current). 
Then substituting ho(x) and j z (x) into the Maxwell equa- 
tions one can obtain the analog of the Ferrell-Prangc 
equation d 2 A9 = \j 2 smA9, where \j = (Kr/2(p ) 1 / 2 
is the Josephson penetration depth. It is convinient to 
introduce new dimensionless coordinate x 1 = X\j2^j nr. 
Then the magnetic part of the Hclmholtz free energy 
F = G + 2BH per unit length can be represented in 
the standard form 

F = G s ^- J |(1 - cosAfl) + i (a^) 2 j dx'. (11) 

From the expression (|11[) it is easy to obtain the value 
of the Josephson critical field h c j which is the minimal 
field of the vortex penetration into the junction^ 3 .. The 
expression for h c j reads as 



(12) 



To describe the temperature dependence of the criti- 
cal magnetic field which corresponds to field penetration 
into the TP we use the results of Ref. HH, where the av- 
eraged magnetization Mj of the Josephson junction and 
the corresponding Hclmholtz energy Fj are calculated as 
implicit functions of the external magnetic field h. Sub- 
stituting these dependencies into the Gibbs free energy 
G of the TP and performing numerical calculations we 
obtain the dependence of the critical magnetic field on 
temperature h c (t) which is shown in Fig. 2 (red solid 
curve). Obviously this dependence exceeds the depen- 
dence hco{t) corresponded to the condition Go = (see 
blue dashed curve in Fig. 2) only for temperatures where 
hco(t) > h c j(t). Note that the obtained dependence h c (t) 
depends only on one fitting parameter r. Thus we hope 
that high-accuracy measurements of the critical magnetic 
fields can allow to estimate the values of the twinning 
plane transparencies for different type I superconductors. 

Note also that the magnetic field penetration into the 
center of the TP can be detected experimentally in the 
Josephson current measurements. Let us consider the 
pair of contacts which are positioned parallel to the TP 
at a distance which is much less than a width of supercon- 
ducting region I cx t -1 / 2 near the TP. Then the Joseph- 
son current through these contacts would be extremely 
sensitive to the magnetic field in the center of the TP. 
Indeed without magnetic field the averaged over the TP 
length Josephson current can be nonzero if AO ^ (see 



(JTOj) ) . Ovcrwise in case of the magnetic field penetra- 
tion the phase difference would increase with the increase 
in x and the corresponding averaged Josephson current 
through the TP would be negligibly small. 

For the type II superconductors the fact of full pen- 
etration of the magnetic field into the center of the TP 
doesn't lead to any substantial consequences since the 
TP in this case weakly screen the external magnetic field. 
To calculate the correction to the magnetic susceptibility 
one can use the approach from Ref. QJ. We will restrict 
ourselves for the case of ultra type II superconductors 
with k > 1, r = oo and the temperatures in the range 
< t < 1. In this case the magnetic field b{z) slightly 
differs from the external field h which allows one to con- 
sider the field profilein the form b(z) = h + 6b(z), where 
|<5&(z)| <C h. Then the value 6b(z) satisfies the equation 
{u 2 - l)d 2 6b = 2h/n 2 , where u = coth {\fi\z\ + p/2) and 
p = ln[(l + \/t)/(l — \/t)}. The boundary conditions in 
the case of full penetration of the magnetic field into the 
TP center are Sb(u = 1) = and 5b(u = i" 1 / 2 ) = 0. 
Note that under the assumption that the magnetic field 
has its minimum on the TP the last condition should be 
replaced with^ d u 5b(u = t^ 1 / 2 ) = 0. The exact solution 
of the equation for 8b allows one to obtain the correction 
to the magnetic susceptibility due to the penetration of 
the magnetic field into the TP. This correction has the 
form (in dimensional units) 



AM, 



(ii) 



1 



K 2 7r(l - Vt) 
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Note that at t — > the correction AM(jn oc ln 2 (t _1 ) 
and is negligibly small since it has weak singularity com- 
pared with the full diamagnetic moment M d ^jj^ of the 
TP, which diverges like t-^Hn 2 ^- 1 ) (see Ref. i). At 
t — >• 1 the correction AM(jj) cx (1 — t) and is also small. 

Thus we have shown that in superconductors with 
twinning planes with low transparency the penetration 
of the parallel magnetic field into the twinning plane is 
energetically favorable. For the type I superconductors 
this leads to the essential increase of the critical magnetic 
field and to the broadening of the temperature range 
where the TPS can exist. Our theory does not contain 
any fitting parameters except r, so this can provide an 
opportunity to estimate the r values for different type 
I superconductors on the basis of the critical magnetic 
field measurements. 
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